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Abstract 


We demonstrate that crystalline color superconductivity may arise as a re¬ 
sult of pairing between massless quarks and quarks with nonzero mass m s . 
Previous analyses of this phase of cold dense quark matter have all utilized 
a chemical potential difference dp to favor crystalline color superconductivity 
over ordinary BCS pairing. In any context in which crystalline color super¬ 
conductivity occurs in nature, however, it will be m s -induced. The effect of 
m s is qualitatively different from that of dp in one crucial respect: m s de¬ 
presses the value of the BCS gap Ao whereas dp leaves Ao unchanged. This 
effect in the BCS phase must be taken into account before m s -induced and 
dp-induced crystalline color superconductivity can sensibly be compared. 
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I. INTRODUCTION 


Since quarks which are antisymmetric in color attract, cold dense quark matter is unsta¬ 
ble to the formation of a condensate of Cooper pairs, making it a color superconductor [jlj . 
At asymptotic densities, the ground state of QCD with quarks of three flavors ( u , d and s ) 
with equal masses is expected to be the color-flavor locked (CFL) phase f|]-£|]. This phase 
features a condensate of Cooper pairs of quarks which includes ud, us, and ds pairs. Quarks 
of all colors and all flavors participate in the pairing, and all excitations with quark quantum 
numbers are gapped. In this phase, left-flavor and right-flavor symmetries are both locked 
to color, breaking chiral symmetry |2J. As in any BCS state, the Cooper pairing in the CFL 
state pairs quarks whose momenta are equal in magnitude and opposite in direction, and 
pairing is strongest between pairs of quarks whose momenta are both near their respective 
Fermi surfaces. Pairing persists even in the face of a stress (such as a chemical potential 
difference or a mass difference) that seeks to push the Fermi surfaces apart, although a stress 
that is too strong will ultimately disrupt Cooper pairing ||||. Thus, the CFL phase persists 
for unequal quark masses, so long as the differences are not too large |5j.[|. This means that 
the CFL phase is the ground state for real QCD, assumed to be in equilibrium with respect 
to the weak interactions, as long as the density is high enough. 

Imagine decreasing the quark number chemical potential fi from asymptotically large 
values. The quark matter at first remains color-flavor locked, although the CFL condensate 
may rotate in flavor space as terms of order m* in the free energy become important [[7]. 
Color-flavor locking is maintained until a transition to a state in which some quarks become 
ungapped. This “unlocking transition”, which must be first order |^J|, occurs when ill! 

h ~ m s/4A 0 . (1) 


In this expression, Ao is the BCS pairing gap, estimated in both models and asymptotic 
analyses to be of order tens to 100 MeV Q, and m s is the strange quark mass parameter. m s 
includes the contribution from any (ss) condensate induced by the nonzero current strange 
quark mass, making it a density-dependent effective mass, decreasing as density increases 
and equaling the current strange quark mass only at asymptotically high densities. At 
densities which may arise at the center of compact stars, corresponding to fi ~ 400 — 500 
MeV, m s is certainly significantly larger than the current quark mass, and its value is not 
well-known. In fact, m s decreases discontinuously at the unlocking transition [|l0] . Thus, 


the criterion (l) can only be used as a rough guide to the location of the unlocking transition 
in nature [[HJ . Given this quantitative uncertainty, there remain two logical possibilities for 
what happens as a function of decreasing fi. One possibility is a first order phase transition 
directly from color-flavor locked quark matter to hadronic matter, as explored in Ref. ||. 
The second possibility is an unlocking transition |3|.[6| to quark matter in which not all quarks 
participate in the dominant pairing, followed only at a lower /i by a transition to hadronic 
matter. We assume the second possibility here, and explore its consequences. 

Once CFL is disrupted, leaving some species of quarks with differing Fermi momenta 
and therefore unable to participate in BCS pairing, it is natural to ask whether there is 
some generalization of the ansatz in which pairing between two species of quarks persists 
even once their Fermi momenta differ. Crystalline color superconductivity is the answer to 
this question. The idea is that it may be favorable for quarks with differing Fermi momenta 
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to form pairs whose momenta are not equal in magnitude and opposite in sign [11,12]. This 


generalization of the pairing ansatz (beyond BCS ansatze in which only quarks with momenta 
which add to zero pair) is favored because it gives rise to a region of phase space where both 
of the quarks in a pair are close to their respective Fermi surfaces, and such pairs can be 
created at low cost in free energy. Condensates of this sort spontaneously break translational 
and rotational invariance, leading to gaps which vary in a crystalline pattern. As a function 
of increasing depth in a compact star, // increases, m s decreases, and Ao changes also. If in 
some shell within the quark matter core of a neutron star (or within a strange quark star) the 
quark number densities are such that crystalline color superconductivity arises, rotational 
vortices may be pinned in this shell, making it a locus for glitch phenomena [[T^ .|Hj|| . 

An analysis of these ideas in the context of the disruption of CFL pairing is complicated 
by the fact that in quark matter in which CFL pairing does not occur, up and down quarks 
may nevertheless continue to pair in the usual BCS fashion. In this 2SC phase, which was 
the earliest color superconducting phase to be studied [I]], the attractive channel involves 
the formation of Cooper pairs which are antisymmetric in both color and flavor, yielding a 
condensate with color (Greek indices) and flavor (Latin indices) structure (q^Qb) ~ 

This condensate leaves five quarks unpaired: up and down quarks of the third color, and 
strange quarks of all three colors. Because the BCS pairing scheme leaves ungapped quarks 
with differing Fermi momenta, crystalline color superconductivity may result. 

To date, crystalline color superconductivity has only been studied in the simplified model 
context with pairing between two quark species whose Fermi momenta are pushed apart by 
turning on a chemical potential difference |[T2[.[T4| [T7| . rather than considering CFL pairing 
in the presence of quark mass differences. Our goal in this paper is to investigate the ways 
in which the response of the system to mass differences is similar to or different from the 
response to chemical potential differences. We can address this question within the two- 
flavor model by generalizing it to describe pairing between massless up quarks and strange 
quarks with mass m s . For completeness, we introduce 


hu = h - 

fjj S — n + 6 /i , (2) 

allowing us to consider the effects of m s and 5/jl simultaneously. We shall use this two-flavor 
toy model throughout, deferring an analysis of crystalline color superconductivity induced 
by the effects of m s on three-flavor quark matter to future work. 


A. Consequences of Sfj, / 0, with m s = 0 


Before describing the consequences of m s ^ 0, let us review the salient facts known about 
the consequences of 5fi ^ o, upon taking rn s = 0. If |<5/i| is nonzero but less than some 5/q, 
the ground state in the two-flavor toy-model is precisely that obtained for 5ji = 0 |T7|JT8|JT^1 .p| 


1 In this two-flavor toy-model the diquark condensate is a flavor singlet. As the condensate breaks 
no flavor symmetries, there is no analogue of the rotations of the condensate in flavor space which 
occur within the CFL phase with nonzero 5/r J7.19]. 
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In this 2SC state, red and green up and strange quarks pair, yielding four quasiparticles with 
superconducting gap A 0 Jl|. Furthermore, the number density of red and green up quarks 
is the same as that of red and green strange quarks. As long as |<5/i| is not too large, this 
BCS state remains unchanged (and favored) because maintaining equal number densities, 
and thus coincident Fermi surfaces, maximizes the pairing and hence the interaction energy. 
As is increased, the BCS state remains the ground state of the system only as long as its 
negative interaction energy offsets the large positive free energy cost associated with forcing 
the Fermi seas to remain coincident. In the weak coupling limit, in which A 0 //i 1, the BCS 

state persists for \S/j\ < 8fii = A 0 /\/2 [|T7].|T^]. For larger A 0 , the 1/^/2 coefficient changes 
in value. These conclusions are the same whether the interaction between quarks is modeled 
as a point-like four-fermion interaction or is approximated by single-gluon exchange [fij. 
The loss of BCS pairing at |<5/i| = 5/ii is the analogue in this toy model of the unlocking 
transition. 

If |h/x| > Sfii, BCS pairing between u and s is not possible. However, in a range Sji\ < 
\8/i\ < 5/j, 2 near the unpairing transition, it is favorable to form a crystalline color supercon¬ 
ducting state in which the Cooper pairs have nonzero momentum. This phenomenon was 
first analyzed by Larkin and Ovchinnikov and Fulde and Ferrell |U] (LOFF) in the context 
of pairing between electrons in which spin-up and spin-down Fermi momenta differ. It has 
proven difficult to find a condensed matter physics system which is well described simply 
as BCS pairing in the presence of a Zeeman effect: any magnetic perturbation that may 
induce a Zeeman effect tends to have much larger effects on the motion of the electrons, as 
in the Meissner effect. The QCD context of interest to us, in which the Fermi momenta 
being split are those of different ffavors rather than of different spins, therefore turns out to 
be the natural arena for the phenomenon first analyzed by LOFF. 

The crystalline color superconducting phase (also called the LOFF phase) has been 
(following Refs. [|ll|]) upon making the simplifying assumption that 


described in Ref. 

quarks interact via a four-fermion interaction with the quantum numbers of single gluon 
exchange. In the LOFF state, each Cooper pair carries momentum 2q with |q| ps 1.25/x. 
The condensate and gap parameter vary in space with wavelength vr/|q|. Although the 
magnitude |q| is determined energetically, the direction q is chosen spontaneously. The 
LOFF state is characterized by a gap parameter A and a diquark condensate, but not by 
an energy gap in the dispersion relation: the quasiparticle dispersion relations vary with the 
direction of the momentum, yielding gaps that vary from zero up to a maximum of A. The 
condensate is dominated by those regions in momentum space in which a quark pair with 
total momentum 2q has both members of the pair within ~ A of their respective Fermi 
surfaces. These regions form circular bands on the two Fermi surfaces. Making the ansatz 
that all Cooper pairs make the same choice of direction q corresponds to choosing a single 
circular band on each Fermi surface. In position space, it corresponds to a condensate which 
varies in space like 


(■0(x)-0(x)) oc Ae 2?qx . (3) 

This ansatz is certainly not the best choice. If a single plane wave is favored, why not two? 
That is, if one choice of q is favored, why not add a second q, with the same |q| but a 
different q? If two are favored, why not three? This question, namely, the determination of 
the favored crystal structure of the crystalline color superconductor phase, is unresolved but 
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is under investigation. Note, however, that if we find a region 8 /ii < \5/i\ < Sfi 2 in which the 
simple LOFF ansatz with a single q is favored over the BCS state and over no pairing, then 
the LOFF state with whatever crystal structure turns out to be optimal must be favored 
in at least this region. Note also that even the single q ansatz, which we use henceforth, 
breaks translational and rotational invariance spontaneously. The resulting phonon has been 
analyzed in Ref. |L5]. 

Crystalline color superconductivity is favored within a window Shi < \S/a\ < 5/i 2 . As 
\Sfj,\ increases from 0, one finds a first order phase transition from the ordinary BCS phase 
to the crystalline color superconducting phase at \S/i\ = 8 Hi and then a second order phase 
transition at |<5/i| = <5/r 2 at which A decreases to zero. Because the condensation energy 
in the LOFF phase is much smaller than that of the BCS condensate at S/i — 0, the value 
of 5Hi is almost identical to that at which the naive unpairing transition from the BCS 
state to the state with no pairing would occur if one ignored the possibility of a LOFF 
phase, namely 5Hi = A 0 /\/2- For all practical purposes, therefore, the LOFF gap equation 
is not required in order to determine 5Hi- The LOFF gap equation is used to determine 
< 5 /i 2 and the properties of the crystalline color superconducting phase |12|. In the limit of 
a weak four-fermion interaction, the crystalline color superconductivity window is bounded 
by 5Hi = Ao/\/2 and <5 /i 2 = 0.754Ao, as first demonstrated in Refs. 0 - These results have 
been extended beyond the weak four-fermion interaction limit in Ref. 0 . 


We now know that the use of the simplified point-like interaction significantly underesti¬ 
mates the width of the LOFF window: assuming instead that quarks interact by exchanging 
medium-modified gluons yields a much larger value of <5/x 2 |L 6 |. This can be understood upon 
noting that quark-quark interaction by gluon exchange is dominated by forward scattering. 
In most scatterings, the angular positions on their respective Fermi surfaces do not change 
much. In the LOFF state, small-angle scattering is advantageous because it cannot scatter 
a pair of quarks out of the region of momentum space in which both members of the pair 
are in their respective circular bands, where pairing is favored. This means that it is natural 
that a forward-scattering dominated interaction like single-gluon exchange is much more 
favorable for crystalline color superconductivity that a point-like interaction, which yields 
s-wave scattering. Thus, although for the present we shall use the point-like interaction in 
our analysis of m s -induced crystalline color superconductivity, it is worth remembering that 
this is very conservative. 


B. Consequences of m s A 0 

In the absence of any interaction, and thus in the absence of pairing, the effect of a 
strange quark mass is to shift the Fermi momenta to 

Pf = h~ Sh 

Pf = + <W - m 2 s . (4) 

Assuming both \Sh\/h an d m s //i are small, the separation between the two Fermi momenta 
is ~ \25h — m 2 s /M- This suggests the conjecture that even when m s 7 ^ 0 the description 
given in the previous subsection continues to be valid upon replacing \Sh\ by |<5/i — m^/4/r|. 
We shall see in Section 3 that this conjecture is incorrect in one key respect: whereas if 
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m s = 0 a 5 /i which is nonzero but smaller than has no effect on the BCS state, the 
BCS gap A 0 decreases with increasing m 2 , as we discuss further in an appendix. We show 
that for small m 2 , A 0 (m s )/A 0 (0) decreases linearly with m 2 . Because A 0 occurs in the free 
energy in a term of order Ag/i 2 , the m s -dependence of A 0 corrects the free energy by of order 
A 0 (0) 2 m 2 . As has no analogous effect, we conclude that m 2 /4p and 5fi have qualitatively 
different effects on the paired state. 

At another level, however, the story is quite similar to that for m s = 0: if \S/j, — m 2 /4p| is 
small enough, we find the BCS state; if \5fi — m 2 /4/i| lies within an intermediate window, we 
find LOFF pairing; if | S/j, — m 2 /4p| is large enough, no pairing is possible. The boundaries 
between the phases, however, are related to a A 0 (m s ), rather than simply to a constant A 0 . 
That is, the definitions of “small enough” and “large enough” are m s -dependent. We map 
the (m s , <5/i) plane in this paper. In Section 2, we derive a gap equation which we then use in 
Section 3 to analyze the BCS phase, finding the region of the plane in which BCS 

pairing is favored, and again use in Section 4 to analyze crystalline color superconductivity, 
finding windows in the (m s ,<5/i) plane in which LOFF pairing is favored. 


II. THE GAP EQUATION 


In this Section, we sketch the derivation of the gap equation which describes either the 
BCS or the LOFF states in a two-flavor theory with m s ^ 0 and <5/i ^ 0. In the crystalline 
color superconducting phase, the condensate contains pairs of u and s quarks with momenta 
such that the total momentum of each Cooper pair is given by 2q, with the direction of 
q chosen spontaneously. As noted above, wherever there is an instability towards (§), we 
expect the true ground state to be a crystalline condensate which varies in space like a sum 
of several such plane waves with the same |q|. (To recover the BCS gap equation from that 
which we now derive, simply set |q| = 0.) 

In order to describe pairing between u quarks with momentum p + q and s quarks with 
momentum — p + q, we use a modified Nambu-Gorkov propagator: 


T(p,g) 


/ ^u(p + q) \ 
ts(p - q) 
ts(~p + q ) 
\^u(-p-q) / 


(5) 


Note that by q we mean the four-vector (0,q). The Cooper pairs have nonzero total mo¬ 
mentum but the ground state condensate fl3|) is static. The momentum dependence of (H) 
is motivated by the fact that in the presence of a crystalline color superconducting con¬ 
densate, anomalous propagation does not only mean picking up or losing two quarks from 
the condensate. It also means picking up or losing momentum 2q. The basis (^]) has been 
chosen so that the inverse fermion propagator in the crystalline color superconducting phase 
is diagonal in p-space and is given by 

fft + i + Hulo 0 -A(p,-q) 0 \ 

0 ]b-fi + ^o-m s 0 A (p,q) . . 

-A (p,-q) 0 (j> - <h - p s 7 o + m s ) T 0 

V 0 A (p,q) 0 (:i> + <k - Pulo) T ) 


S \p,q ) = 
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where A = yoA^yo and the matrix A is given by A = AC^e 01 ^ 3 . Note that the condensate 
is explicitly antisymmetric in flavor. 2p is the relative momentum of the quarks in a given 
pair and is different for different pairs. In the gap equation below, we shall integrate over p 0 
and p. As desired, the off-diagonal blocks describe anomalous propagation in the presence 
of a condensate of diquarks with momentum 2q. 

We obtain the gap equation by solving the one-loop Schwinger-Dyson equation given by 

S~\k,q) - So\k,q) = -g 2 j -^-V A S{p,q)V^ D% B (k - p), (7) 


where D^ B = D ,ju 8ar is the gluon propagator, S is the full quark propagator, whose inverse 
is given by (H), and So is the fermion propagator in the absence of interaction, given by S 
with A = 0. The vertex is defined as 

/ 7/A a /2 0 0 

_ 0 7m A a /2 0 

11 o 0 -( 7/ A a /2) t 

V o o o 


0 
0 
0 

-^\ A m T ) 


( 8 ) 


As in Refs. [ |H2| , |T4| , we simplify the interaction between quarks by making it point-like. 
That is, we replace g 2 D^ v by 3 Gg^. Our model Hamiltonian has two parameters, the 
four-fermion coupling G and a cutoff A on the three-momentum p. The dependence on 
the cutoff has been studied in Ref. [ 121: upon varying A while at the same time varying 


the coupling G such that Ao is kept fixed, the relation between other physical quantities 
and Ao is reasonably insensitive to variation of A. We often use the value of the BCS gap 
A 0 obtained at m s — Sp — 0 to parametrize the strength of the interaction. We can do 
this because A 0 increases monotonically with increasing G. Doing so proves convenient 
because both 5p\ and 5p2 are given simply when written in terms of the physical quantity 
A 0 , whereas writing them in terms of the model-dependent parameters G and A requires 
unwieldy expressions. 

After some algebra (essentially the determination of S given S~ 3 specified above), and 
upon suitable projection, the Schwinger-Dyson equation (|7|) reduces to a gap equation for 
the gap parameter A given (in Euclidean space) by 


A 


d 4 p _ 8GA(A 2 + w 2 ) _ 

(27t) 4 A 4 + 2A 2 w 2 + [(po + ip + i5p) 2 — |p — q| 2 — m 2 } [(p 0 — ip + iSp) 2 


p + ql 2 ] 
( 9 ) 


where 


2 _ 2 i 2 

W = p + p 


9 2 + (Po + i<V) ! 


( 10 ) 


Eq. (j9|) reduces to the LOFF gap equation of Ref. [[L4|] when m s = 0. The gap equation ([|) 


includes the contributions of pairing involving antiparticles, obviously far from the Fermi 
surfaces. These contributions can be removed using the methods of Ref. |14| (and if one 
does so and then sets m s = 0 one obtains the gap equation of Ref. |P~2|). However, as 
the antiparticle contributions are suppressed by A /p and as removing them leads to rather 
unwieldy expressions, we shall simply use eq. (f|). 
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III. THE BCS GAP AND THE BCS REGION 


We expect BCS pairing to be possible when the u and s Fermi momenta in noninteracting 
quark matter are similar in value. From (|4|) we see that p F = p s F in noninteracting quark 
matter when 

, m 2 * . . 

5p = — . 11 

4/i 

Note that we have not assumed that dp or m s is small compared to p in this expression. 
Let us now analyze the BCS gap equation, keeping this in mind. Upon setting q = 0, Eq. 
(|) reduces to the BCS gap equation of Ref. ||: 

A = f d 4 p _ 8GA 0 (A 2 + w 2 ) _ 

0 J (27t) 4 Aq + 2Aq Wq + [(po + ip + idp) 2 - p 2 - m 2 s ] [(p 0 - ip + idp) 2 - p 2 } 


where 


w? 


= p 2 + p 2 + (po + idpf 


and where we use the notation A 0 for the BCS gap, to keep track of the fact that here q = 0. 
Eq. (|T^) involves an integration over p 0 from — oo to +oo, which can be thought of as a 
contour in the complex p 0 -pl ane ; which we choose to close in the upper half-plane. With 
q = 0, the poles in the integrand of Eq. (jT2|) are easily analyzed. There are four poles in 
the po _ pl ane - Of these, two correspond to the contribution of antiparticle pairing and have 
residues which are suppressed by A 0 //i relative to the residues of the two remaining poles, 
which correspond to particle pairing. The residues of the two particle poles are equal in 
magnitude and opposite in sign. As a result, the existence of a nonzero Ao which solves ( 12 ) 
requires that one and only one of the particle poles lies in the upper half-plane. If we choose 
dp according to 0 . we do in fact End one particle pole in the upper half-plane and thus a 
nonzero, m s -dependent, solution A 0 as shown in Fig. [I]. A 0 decreases with m 2 , linearly for 
small m 2 . The m s -dependence of A 0 can be further understood, but as this would divert us 
we defer doing so to an appendix. 

What is the effect of changing dp away from the choice (|ll])? This is simply understood 
upon noticing that po and dp enter the gap equation only in the combination p 0 + idp 
This means that the only effect of dp is to shift the location of the poles. Whereas m s 
affects the value of the residues, dp does not. As long as |dp — m 2 /Ap\ < A o ,0 there is one 
particle pole in the upper half-plane and A 0 is identical to that found when dp = m 2 s /Ap, as 
in Fig. p]. Once |dp — m 2 /4:p\ > A 0 , there are either zero or two particle poles in the upper 
half-plane and therefore no such solutions A 0 exist. 

As long as |dp — m 2 /Ap\ < A 0 (m s ), we have found a dp- independent, m s -dependent 
solution Ao (m s ) (shown in Fig. |T|) that describes a BCS state which is a local minimum of 
the free energy 0 at a given m s and dp. However, this BCS state need not be global minimum 


2 The criterion is \5p — m 2 /4p\ < Ao in the limit in which all these quantities are small compared to 

p. As we shall see momentarily, this is not the most important criterion anyway, and we therefore 

do not specify it precisely. 
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FIG. 1. BCS gap Ao vs. m^/4/i, with 5fi set equal to m 2 / 4/i to ensure that in the absence of 
any pairing, the u and s Fermi momenta would be equal. In this and in all subsequent figures, we 
set /i = 400 MeV and choose a cutoff A = 1000 MeV. We have chosen three different values of the 
coupling constant G, yielding Ao = 10,40,80 MeV at m s = 0. Ao decreases with increasing m 2 , 
linearly for small m 2 . 


of the free energy. To determine the region in the (mg/4/i,5/i) plane in which BCS pairing 
is favored over impaired quark matter, we must evaluate the free energy difference between 
the BCS state and the unpaired state. This is given by integrating the gap equation over A 
up to the equation’s solution, A 0 (m s ). 


n 


BCS 


n 


unpaired 


rA 0 (m a ) / A 

L G 


AF integrand) , 


(13) 


where “integrand” refers to that of the right-hand side of Eq. ([12|) . The curves in the 
(m^/4/i, 5/j.) plane along which the above expression vanishes determine the region in which 
the BCS phase is favored. These curves are plotted in Fig. ^]. At any given rn s , BCS pairing 
is favored within a range of values of around m 2 / 4/i. Everywhere within this range, 
A 0 takes on the same (m s -dependent; 5/i-independent) value. The boundaries of the BCS 
region are given approximately by 


, rn; 

4/i 


< 


A 0 (m s ) 

V2 


(14) 


derived in Ref. [|J and valid if A 0 (m s ), rn s and 5fi are all small compared to /i. We see that 
as m s increases and A 0 (m s ) decreases, the BCS region narrows. 










100 


80 

60 

Sp- 40 

(MeV) 

20 


-20 

m l / 4 p (MeV) 

FIG. 2. The region in the (m^/A/i, <5p) plane where the BCS state is favored over unpaired 
quark matter. We have chosen a single value of the coupling, corresponding to Ao(0) = 40 MeV. 
At each value of m^/4/q Ao (m s ) is independent of 5fi within the shaded region and zero outside 
it. Ao(m s ) decreases with m 2 /4p as shown in Fig. |l]. 



IV. THE CRYSTALLINE COLOR SUPERCONDUCTIVITY WINDOWS 


We now determine the regions of the (m;?/4/i, 5/j.) plane where the crystalline color su¬ 
perconducting phase, with |q| 7 ^ 0, is favored. These regions will be strips just above and 
just below the BCS region of Fig. |j. Let us denote the upper boundary of the BCS region 
5fii(m s ). We now show that at a given m s , the crystalline color superconducting state is 
favored in some region 5^i{rn s ) < 5ji < dfi 2 (m s )fl There is also a crystalline color super¬ 
conductivity window of m s -dependent width just below the BCS region. 

At 5/i = h/i 2 , the gap A in the LOFF state goes continuously to zero. This means 
that we can determine h/i 2 as follows. We set A = 0 in the gap equation ([|), obtaining a 
“zero-gap relation” among <5/i, rn s , and |q|. We choose a value of m s , and find the largest 
value of 5/j, at which the zero-gap relation can be satisfied. This is h/i 2 . (The smallest 
value of h/i at which the zero-gap relation can be satisfied is the lower boundary of the 
crystalline color superconductivity window below the BCS region of Fig. 1) In Fig. | 
we plot [h/i 2 (m s ) — 5/ii(m s )]/A 0 (m s ) vs. m 2 / 4/i for four different values of the coupling, 
corresponding to A 0 (0) = 10,40,80,100 MeV. At m s = 0 and in the weak coupling limit 
in which A 0 (0) —■> 0, [<5/t 2 (0) — h/q(0)]/A 0 (0) —■> [0.754 — 0.707] = 0.047 [Jill. The results 
at m s = 0 with A 0 7 ^ 0 which we find in Fig. p] agree with those derived in Ref. |T 2 


These results are symmetric about h/i = 0. For m s 7 ^ 0, we now see that when it is scaled 
by A 0 (m s ), the width of the crystalline color superconductivity window changes little, for 


3 The value of <5/ri(?n s ) was determined in the previous section by comparing the free energy of 
the BCS and unpaired states. Here, we should in principle compare BCS and LOFF states to 
determine 5pi(m s ), but this makes very little difference jf^]. 


9 












nr/ 4 p (MeV) 


FIG. 3. Upper panel: width of the crystalline color superconductivity window above the BCS 
region. At a given m s , crystalline color superconductivity occurs for 5[i\{rn s ) < 5(jl < 5^2 (m. s ). 
Here, we show the width of the window 8 / 12 (m s ) — relative to Aq(to, 5 ) for four val¬ 

ues of the coupling. The four curves from top to bottom correspond to couplings for which 
A 0 (0) = 10,40,80,100 MeV. Lower panel: width of the crystalline color superconductivity window 
below the BCS region. The curves from top to bottom correspond to increasing couplings, as in 
the upper panel. 

smaller values of Ao(0). For larger values of Ao(0), the width of the window increases as m s 
increases from zero. Indeed, we see that for Ao(0) = 100 MeV, the window has closed at 
m s = 0 whereas crystalline color superconductivity continues to occur at large enough m s . 

What do we conclude from our results? First, it is clear that the appropriate variable to 
use to describe the width of the crystalline color superconductivity window is Sfj,/A 0 (m s ), 
as opposed to <5/i/A 0 (0). When plotted as in Fig. [J[. the width of the window is almost 
independent of m s at weak coupling, indicating that we have found the correct scaling 
variable. This means that if we were to plot the crystalline color superconducting windows 
as strips above and below the BCS region of Fig. || the “horizontal width” of the window 
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on the dp — 0 axis, expressed as a window in m^/[4/iA 0 (m s )], is the same as the “vertical 
width” of the windows in dp/A 0 (0) on the rn s = 0 axis. Thus, to the extent that the 
curves in Fig. ^ are horizontal we conclude that m, s -induced and dp-induced crystalline color 
superconductivity are equally robust. However, the curves in Fig. ^jare not quite horizontal. 
At all but the weakest of couplings, the width of the crystalline color superconductivity 
window increases with m s , meaning that crystalline color superconductivity is somewhat 
more robust if it is m, s -induced than if it is dp-induced. Indeed, we find that at the moderate 
coupling corresponding to A 0 (0) = 100 MeV, m s -induced crystalline color superconductivity 
occurs whereas dp-induced crystalline color superconductivity does not. 
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m ^ / 4 p(m s ) (MeV) 



FIG. 4. BCS gap Ao vs. m 2 /4p(m s ), with p(m s ) taken to increase with increasing m s in 
just such a way that, in the absence of interaction, the Fermi momenta of both the u and s 
quarks are m s -independent. We plot the ratio Ao(m s )/Ao(0), for couplings corresponding to 
Ao(0) = 10,40,80 MeV, from bottom to top. The curves are essentially independent of the 
coupling when plotted in this way. 

APPENDIX A: THE M^-DEPENDENCE OF A 0 

In this appendix, we seek to better understand the m s -dependence of the BCS gap A 0 , 
shown in Fig. [lj. To the extent that the coupling can be taken as weak when Ao = 80 MeV, 
one might expect that when each of these three curves are scaled by their value of Ao(0), 
they should fall on top of one another. This does not occur. The reason for this can be 
understood upon first setting the interaction to zero. In the absence of pairing, increasing 
rn s while setting Sp = m 2 / 4/i to ensure that p F = p s F has two effects: the value of pp 
decreases, and the strange quark Fermi velocity decreases. Since two different effects are 
at play, it is not surprising that the curves in Fig. [L] do not scale trivially. To disentangle 
these effects, we have redone the analysis of Fig. |T] as follows. As we increase m s , we now 
increase p , while keeping Sp = m 2 /Ap(m s ). As before, the latter condition ensures that, in 
the absence of interaction, p F = p s F . Here, we choose p{m s ) in just such a way that both p F 
and p s F are m s -independent. Any remaining m s -dependence of A 0 can now be seen as arising 
from the variation of the difference between the Fermi velocities in the unpaired state. The 
plot analogous to Fig. [I] is given in Fig. £|. and the scaling is now manifest: the curves are 
almost independent of coupling. Note also that A 0 decreases linearly with m 2 as before, but 
the slope in Fig. [4] is much less than that in Fig. [T]. This indicates that the m. s -dependence 
in Fig. |l] can be attributed more to the variation of p F in the unpaired state than to the 
variation of v F in the unpaired state. 
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